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ABSTRACT 



The middle-aged PSR J0357+3205 is a nearby, radio-quiet, bright 7-ray pulsar 
discovered by the Fermi mission. Our previous Chandra observation revealed a 
huge, very peculiar structure of diffuse X-ray emission, originating at the pulsar 
position and extending for > 9' on the plane of the sky. To better understand the 
nature of such a nebula, we have studied the proper motion of the parent pulsar. 
We performed relative astrometry on Chandra images of the field spanning a 
time baseline of 2.2 yr, unveiling a significant angular displacement of the pulsar 
counterpart, corresponding to a proper motion of (X'165 ± 0'.'030 yr -1 . At a 
distance of ~ 500 pc, the space velocity of the pulsar would be of ~ 390 km s" 1 
assuming no inclination with respect to the plane of the sky. The direction of the 
pulsar proper motion is perfectly aligned with the main axis of the X-ray nebula, 
pointing to a physical, yet elusive link between the nebula and the pulsar space 
velocity. No optical emission in the Hq, line is seen in a deep image collected at 
the Gemini telescope, which implies that the interstellar medium into which the 
pulsar is moving is fully ionized. 

Subject headings: Stars: neutron — pulsars: general — pulsars: individual (PSR 
J0357+3205) 
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Introduction 



The Large Area Telescope onboard the Fermi satellite (Fermi-L AT , lAtwood et al 



20091 ) has opened a new era for pulsar astronomy, by detecting 7— ray pulsations (at 



E>100 MeV) from more than 120 pulsars^, about 30% of which are not detected at radio 
wavelengths. The middle-aged PSR J0357+3205 (characteristic age tq ~ 0.54 Myr) is 
one of the most i nteresting radio-qui et pulsars discovered in blind periodicity searches in 



Fermi-LAT data 
bright source list. 



Abdo et al. 



Abdo et al 



2009aj). Its high 7-ray flux (it is included in the Fermi-LAT 
2009bl ). low spin-down luminosity (E rot = 5 x 10 33 erg s -1 ) and 
off-plane position (Galactic latitude b ~ —16°) point to a small distance of about 500 pc. 



We investigated the field of PSR J0357+3205 with a joint X-ray and optical program with 
Chandra and the NOAO Mayall 4m telescope at Kitt Peak. This allowed us to identify the 



soft X-ray counterpart of the pulsar as an u nremarkable source, 



than other well known middle-aged pulsars (IDe Luca et al. 



ooking fainter (and colder) 



201lh . More interestingly, our 



deep Chandra observation unveiled the existence of a very large, elongated feature of diffuse 
X-ray emission, apparently originating at the pulsar position and extending for more than 9' 
(corresponding to ~ 1.3 pc at the distance of 500 pc, assuming no inclination with respect 
to the plane of the sky), with a hard spectrum consistent both with a power law and with a 
hot thermal bremsstrahlung. 



Elongated "tails " of diffuse emission 



nave 



powered puls ars (e.g. 



nebulae (see 



Kargaltsev fc 



Gaensler fe Slane 



'avlov 



Deen associated with several rotation- 



20081 ) and explained as bow-shock pulsar wind 



20041 for a review), where their elongated morphology 
is "velocity-driven". Indeed, if the pulsar moves supersonically through the interstellar 
medium, the termination shock of the pulsar wind assumes a "bullet" morphology, due to 
ram pressure. Particles accelerated at the shock emit synchrotron radiation and cool down, 



1 See https : //confluence . slac . Stanford . edu/display/GLAMCOG/Public+List+of +LAT-Detected+< 
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confined by ram pressure in an elongated region aligned with the pulsar space velocity. 



However, explaining the nature of the nebul a associated wi t 
out to be very challenging. As discussed by 



h PSR J0357+3205 turned 



De Luca et al 



( 1201 ll ). the standard picture 



cannot apply here since the morphology is very different from the "cometary" shape which 
characterizes all other X-ray bow-shock nebulae. There is no emission in the surroundings 
of the pulsar, where the brightest portion (the termination shock) should be - indeed, the 
surface brightness grows as a function of the distance from PSR J0357+3205. Moreover, 
there are no evidences for spectral evolution as a function of the position, at odds with 
expectations for a population of particles injected at the shock and cooling via synchrotron 
radiation. 

Other pictures could be explored. For instance, PSR B 2224+65, the fast moving; 



pulsar powering the well known "Guitar nebula" seen in H Q (ICordes et al. 



19931 ). displays 



an elongated X-ray feature which is reminiscent of the one of our target and cannot be 
a bow-shock nebula because it is misali gned by ~ 118° with respect to the direction of 
the proper motion (IHui fc Becker! 120071 ). Thus, the possibility of a ballistic jet (similar 
to Active Galactic Nuclei), or the hypothesis of a nebula confined by a pre-existing, 



large scale magnetic fie 



Johnson fc Wang 



2010 



d in the interst ellar medium have been proposed (IBandiera 



HuietaL 



2012I ). 



2008 



Indeed, a crucial piece of information in order to understand the physics of the huge 
elongated feature associated with PSR J0357+3205 is the direction of the pulsar proper 
motion. Detecting a pulsar angular displacement aligned with the nebula's main axis would 
link the morphology of the diffuse structure to the pulsar velocity. Conversely, if it were 
misaligned, the case of PSR J0357+3205 would become very similar to the one of PSR 
B2224+65 and would require a different explanation for the nature of the nebula. 



Usually, a pulsar proper motion is measured in the radio band, or, more rarely, in the 
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optical domani. Unfortunately, our target is radio-quiet and has no optical counterpart; 
moreover, timing analysis of gamma-ray photons is not particularly sensitive to the proper 
motion (position al accuracy based on 5 yr of Fermi-LAT timing is estimated to be ~ 2", 



Ray et al. 



20111 ) . The only way to search for a possible proper motion rests on the 
comparison of multi-epoch, high-resolution X-ray images. To this aim, we have obtained 
a multi-cycle observing campaign with Chandra, consisting of two observations to be 
performed at the end of 2011 and at the end of 2013. We will report here on the first 
observation of our program, as well as on a very recent observation of the field in the H Q 
band performed with GMOS instrument at the Gemini North telescope. Indeed, pulsars 
moving supersonically into warm interstellar gas can generate optical emission in the H a 
line, due to collisional excitation of neut ral hydrogen and c harge exchange occurring at 



(and behind) the forward shock (see e.g. 



Cordes et al. 



19931 ). yielding a limb-brightened, 



arc-shaped bow-shock nebula, located at the apex of the forward shock, in the direction of 
the proper motion. 



2. Measurement of the pulsar proper motion 

The superb angular resolution of the Chandra optics makes it possible to measure tiny 
angular displacements of an X-ray source by performing relative astrometry on multi-epoch 
images. Indeed, such an approach has already been used to measure the proper motion 



(De Luca et al. 


2009 


, as well as t 


Kaplan et al. 


C 


20091), 


Becker et al. 



(2012). 



Motch et a 



van Etten et al. 



(2007) 



fl2012f ). 



Motch et al. 



(120081 ) 
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2.1. Chandra observations and data reduction 

Our new observation of PSR J0357+3205 with Chandra was performed on 2011, 
December 24 (Obs.Id. 14007, 29.4 ks exposure time - hereafter tagged as "2011"). Previous 
observations were performed on 2009 October 25 (obs. id. 12008, 29.5 ks - hereafter 
"2009a") and on 2009 October 26 (obs. id. 11239, 47.1 ks - hereafter "2009b"). All data 
were collected using the Advanced CCD Imaging Spectrometer (ACIS) instrument in Timed 
Exposure mode with the VFAINT telemetry mode. We retrieved "Level 1" data from the 
Chandra Science Archive and reprocessed them with the chandra-repro^ script of the 
Chandra Interactive Analysis of Observation Software (CIAO v4.4)[]. 

For each observation, we generated an image in the 0.3-8 keV energy range using 
the original ACIS pixel size (O'.'492/pixel). We performed a source detection using the 
wavdetecto task with wavelet scales ranging from 1 to 16 pixels with a \J2 step size, 
setting a detection threshold of 10 -6 . In all observations, the target was imaged close to the 
aimpoint, on the backside-illuminated chip S3 of the ACIS-S array. 

We cross-correlated the resulting source lists using a correlation radius of 3" and 
we extracted a catalogue of common sources for each pair of observations. In view of 
the density of sources in each image, the possibility of a chance alignement of two false 
detections is < 10~ 5 . As a further step, we selected sources within 4' of the aimpoint since 
the telescope point spread function deteriora tes as a function of offaxis angle, hampering 



source localization accuracy (see discussion in 



De Luca et al 



20091 ). Such an exercise yielded 



12 common sources (2011 vs. 2009a), 10 sources (2011 vs. 2009b) and 16 sources (2009a vs. 



2 http : / / cxc . harvard . edu /ciao/ ahelp / chandr a_repro . html 



3 http : / / cxc . harvard . edu /ciao/ index, html 



http: / / cxc.harvard.edu/ciao/threads/wavdetect / 
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2009b), in addition to the target pulsar. The uncertainty on the source localization on each 
image depends on the signal to noise as well as on the off axis angle and ranges from ~ 0.08 
to ~ 0.7 pixel per coordinate. 



2.2. Relative astrometry 

The positions of the selected, common sources (excluding the pulsar counterpart) were 
adopted as a reference grid to perform relative astrometry. We used the ACIS SKY reference 
systeno (pixel coordinates with axes aligned along Right Ascension and Declination). 
Taking into account the corresponding uncertainties, we computed the best geometric 
transformation needed to superimpose the reference frames of two images collected at 
different epochs. 

We superimposed the most recent 2011 data to first-epoch 2009a and (separately) 
2009b data in order to measure the possible pulsar displacement over a baseline of ~ 2.2 
yr. We also superimposed 2009a data to 2009b data (no displacement is expected on a 
1-day baseline) in order to check for systematics affecting our analysis. A simple translation 
yielded a good superposition in all cases (see Table 1). The uncertainty on the frame 
registration turned out to be smaller than 50 mas per coordinate. Adding a further 
free parameter to the transformation (a rotation angle) did not result in a statistically 
compelling improvement of the fit. 

We applied the best-fit transformation to the coordinates of the pulsar counterpart and 
we computed its displacement between different epochs. The error budget for the overall 
pulsar displacement includes the uncertainty on the pulsar position in each image as well 
as the uncertainty on the multi-epoch frame registration. Results are shown in Table 1. 



5 http:/ /cxc. harvard.edu /ciao/ ahelp / coords, html 
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Table 1: Results of X-ray relative astrometry 





2011 vs. 2009a 


2011 vs. 2009b 


2009a vs. 2009b 


Time baseline 


2.16 yr 


2.16 yr 


1 day 


Number of ref.srcs 




10 


16 


uncertainty on X shi f t (pixels) 
X 2 (dof) 


0.09 
13.6 (10) 


0.08 
15.6 (9) 


0.06 
7.8 (15) 


uncertainty on Y s hift (pixels) 
X 2 (d.o.f.) 


0.08 
8.1 (10) 


0.07 
13.0 (9) 


0.06 
13.8 (15) 


PSR X displacement (pixels) 
PSR Y displacement (pixels) 


0.53 ±0.12 
0.54 ±0.11 


0.50 ±0.11 
0.50 ±0.10 


0.10 ±0.10 
0.04 ±0.10 







a A reference source yielding high residuals was rejected. 
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Displacement of the pulsar as well as residuals on the positions of the reference sources after 
frame registration are also shown in Figure 1. A significant and consistent displacement of 
the pulsar is apparent both in the 2011 vs. 2009a and in the 2011 vs. 2009b comparison, 
while no displacement is seen for the pulsar in the 2009a vs. 2009b one. 

By combining the two measurements, we obtain a pulsar proper motion fi a cos(5) = 
117 ± 20 mas yr -1 and fis = 115 ± 20 mas yr _1 . This corresponds to a total proper motion 
of 165±30 mas yr" 1 , translating to a (projected) space velocity of ~ 390d 50 o km s" 1 (where 
dsoo is the distance to the pulsar in units of 500 pc). The position angle of the proper 
motion is 314° ± 10° (North to East) (see Figure 2). The proper motion direction is very 
well aligned with the main axis of the nebula, which is indeed an "X-ray trail" . 

3. Gemini H Q observations. 

We observed the field of PSR J0357+3205 in the H a band on September 23, 2012 
with the Gemini-North telescope on the Mauna Kea Observatory. We used the Gemini 
Multi-Object Spectrograph (GMOS) in its imaging mode equipped with the interim upgrade 
e2v deep depletion (DD) detector, with enhanced response in the blue and red arm of the 
spectrum with respect to the original EEV detector. The e2v DD detector is an array of 
three chips, with an unbinned pixel size of 0'.'072 8 and covers an unvignetted field-of-view 
of 5(5 x 5(5. We used the Ha_G0310 filter (A = 656 nm; AA = 7 nm), centered on the H a 
rest wavelength. The e2v DD detector was set in the 2x2 binning mode and read through 
the standard slow read-out/low gain mode with the six amplifiers. Observations were 
performed with an average airmass of 1.24, image quality of ~ 0'/8, and grey time. Two sets 
of six 500 s exposures were obtained, for a total integration time of 6000 s. Exposures were 
dithered in steps of ± 5 pixels to achieve a better signal-to-noise (S/N). 
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We processed the GMOS images using the dedicated CMOS image reduction package 
available in IRAF. After downloading the closest bias and sky flat-field frames from the 
Gemini science archived, we used the tasks gbias and gif lat to process and combine the 
bias and flat-field frames, respectively. We then reduced the single science frames using 
the task gireduce for bias subtraction, overscan correction, image trimming and flat-field 
normalisation. From the reduced science images, we produced a mosaic of the three GMOS 
chips using the task gmosaic and we used the task imcoadd to average-stack the reduced 
image mosaics and filter out cosmic ray hits. We computed the astrometry calibration using 
the wcstoo/J^ suite of programs, matching the sky coordinates of stars selected from the 



Two Micron All Sky Survey (2MASS) All-Sky Catalog of Po int Sources (ISkrutskie et al. 



1996|). 



20061 ) with their pixel coordinates computed by Sextractor ( jBertin fc Arnouts 
After iterating the matching process applying a cr-clipping selection to filter out obvious 
mismatches, high-proper motion stars, and false detections, we obtained mean residuals of 
~ 0'.'2 in the radial direction, using up to 30 bright, but non-saturated, 2MASS stars. To 
this value we added in quadrature the uncertainty a tr = 0'.'08 of t he image registration on 



Lattanzi et al. 



19971 ) 



the 2MASS reference frame. This is given by a tr =^/n/Nsas (e.g. 
where N$ is the number of stars used to compute the astrometric solution, n=5 is the 
number of free parameters in the sky-to -image transformatio n model, as ~ 0'.'2 is the 



mean absolute position error of 2MASS (ISkrutskie et al 



20061 ) for stars in the magnitude 



range 15.5 < K < 13. After accounting for t he O'/OIS uncertainty on the link of 2MASS to 



the International Celestial Reference Frame (ISkrutskie et al. 



20061 ). we evaluated with the 



overall accuracy on our absolute astrometry to be of ~0'/22. 



Unfortunately, no observations of spectro-photometric standards stars were available 



6 http:/ /cadcwww.dao.nrc.ca/gsa/ 



7 http://td c-www.harvard.edu/wcstools/ 
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for the flux calibration of our H a image. Thus, we cross-correlated instrumental magnitudes 
of more than 150 non-saturated stars in the Gemini image to their Rf magnitudes listed in 



the Guide Star Catalogue v2.3.2 (GSC2, 



Lasker et al 



20081 ). which yielded a rather good fit 



with a r.m.s. of 0.16 mag. To assess the GSC2 flux zeropoint, we cross-correlated the GSC2 
Rf magnitude of more than 500 stars with their R magnitude as tabulated in the Stetson 
Standard photometric star archivq£|. We evaluated the transformation as _Rp=0.97R-0.32 
(x 2 = 617, 520 d.o.f.). Then, we computed the H Q fluxes of the 150 selected stars on our 
Gemini image using specific fluxes corresponding to their Rp magnitudes and assuming a 
flat spectrum within the R filter bandpass. This yielded a flux calibration of the H Q image 
with an uncertainty of ~ 0.2 mag. 

No structures of diffuse emission unambiguously related to the fast motion of PSR 
J0357+3205 can be discerned on our Gemini image - in particular, no arc-shaped 
"bow-shock" nebula is seen at the expected position in the direction of the pulsar proper 
motion (see Figure 3). We measured the background properties in the region surrounding 
the pulsar position. Considering a region of 10 square arcsec as a reference and taking into 
account the uncertainty in the flux calibration of our H a image, we computed the 5a upper 
limit to the surface brightness of an unseen bow-shock nebula to be of 5 x 10~ 18 erg cm -2 
s _1 arcsec -2 . 



4. Discussion 

Relative astrometry on our multi-epoch Chandra images unveiled a significant proper 
motion of 165 ± 30 mas yr _1 for PSR J0357+3205, corresponding to a velocity of ~ 390dsoo 
km s _1 , along a direction almost coincident with the main axis of the elongated X-ray 

8 http : //www3 . cadc-ccda.hia-iha.nrc-cnrc .gc . ca/community/STETSON/ar chive/ 
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nebula. 



Our result points to a direct link between the nebula morphology (and physics) and 
the space velocity of the pulsar - supersonic for any reasonable condition of the interstellar 
medium (typical values are of ~ 1, ~ 10 and ~ 100 km s _1 for the cold, warm and 



hot components, respectively, see e.g. 



Kulkarni fc Heiles 



19881 ). As already mentioned, 



"velocity-driven" pulsar wind nebulae are a well-known astrophysical reality. Could the 
X-ray trail of PSR J 0357+3205 be explained as a bow-shock PWN? As discussed by 
De Luca et al.l ( 12011 ). the lack of any evidence of the pulsar wind termination shock could 
be due to its confusion with the pulsar emission. Assuming a velocity of 390d5oo km s _1 and 
no inclination with respect to the plane of the sky, the expected angular separation between 



the pulsar a nd the forward and 



respectively ( IDe Luca et al. 



backward termination shock is $!3n 1 ISM and < 2"n 1ISM , 
1201 ll ). where uijsm is the density of the interstellar medium in 
units of 1 particle cm -1 . Thus, high inclination of the pulsar velocity and/or high density of 
the ISM, and/or a distance larger than expected should be assumed; any of these hypotheses 
could be possible (though rather unlikely). In any case, other peculiarities of the trail 
phenomenology cannot be accounted for in the standard scenario of synchrotron emission 
from shocked pulsar wind. One should explain why the maximum luminosity is generated 
at a distance as large as ~ 3 light years from the parent neutron star. This would require 
an increase of the density of the radiating particles on the same distance scale (possibly due 
to deceleration of the bulk flow?), or an increase of the magnetic field intensity, or a change 
of the angle between the magnetic field and the particle flow. Moreover, the lack of any 



-1/2 



spectral steepening across the trail (IDe Luca et al. 



20 111 ) clashes with the expected, severe 



synchrotron energy losses of particles radiating in the keV range. Thus, either some very 
peculiar phenomenon is occurring in the magnetized pulsar particle outflow (local particle 
re-acceleration? large-scale organization of the flow magnetic field?) or the nature of th e 



nebula is different (but linked to the pulsar velocity). For instance, 



Marelli et al 



(120121 ) 



13 



propose that the trail could be due to thermal emission by the ISM heated by the shock 
driven by the fast moving pulsar. 

The measure of the pulsar space velocity can be used, together with the upper limit 
to the surface brightness of any undetected diffuse structure in our Gemini H Q image, 
to constrain the conditions of the medium in which PSR J0357+3205 is moving. The 
expected flux of the narrow-line component of a pulsar H a bow-shock nebula depends 
on the pulsar spin-down luminosity E rot , the pulsar space velocity v psr , the distance to 
the source Dy 9ri and the neutra l fra ction of the medium X^m- Using the scaling law 



proposed by 



Cordesetal 



( 11993( 1 and 



Chatterj ee Corded (120021 ). we estimated that for any 



reasonable value of the distance to PSR J0357+3205 (ranging from ~ 200 pc, as suggested 
by the non negligible X-ray absorbing column, to ~ 900 pc, where 7-ray luminosity would 
exceed the spin-down luminosity) and in a broad range of possible space velocities for the 
pulsar (assumining the inclination angle with respect to the plane of the sky in the ±75° 
range), the lack of detection of a Ha nebula can only be ascribed to a neutral fraction 
X/sAf <0.01, i.e. the gas surrounding the pulsar is fully ionized. Thus, PSR J0357+3205 
could be traveling across a region filled with ISM in the hot phase. Some contribution to 
the ionization of the medium from the pulsar itself is also p ossible, although its thermal 



emission is not particularly prominent ( iDe Luca et al. 



201lh . 



PSR J0357+3205 is located well away from the Galactic plane (b~ —16°). In view 
of this, it is interesting to note that the direction of the pulsar proper motion is almost 
parallel to the Galac tic plane (it has an inclination of ~ 2° towards the Plane). As done in 



Mignani et al. 



(120121 ). we used the proper motion information to extrapolate back in time 



the trajectory of the pulsar in the Galactic potential to find 



possible associations with open 



clusters and OB associations (iDias et al 



20021 : 



Melhik et al. 



2009 



de Zeeuw et al. 



1999j), 



leaving the pulsar radial velocity as a free parameter and assuming a distance of 500 pc 
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and an age of 0.54 Myr. However, we found no association with clusters or OB associations 
closer than ~ 2 kpc. In any case, such an exercise showed us that if PSR J0357+3205 has a 
radial velocity of order 600 km s _1 (or greater) in the direction away from the Solar System, 
the pulsar could have been born on the Galactic plan e (assuming a scale height of 60 pc, 



typical for massive stars, see e.g. 



Maiz-Apellaniz 



200ll ). Thus, the h ypothesis of a 



r unawa y 



high-mass star as the progenitor of PSR J0357+3205, suggested by iDe Luca et al. 
is not required in a range of possible orbital solutions. 



(2011), 



5. Conclusions 

The radio-quiet PSR J0357+3205 is traveling very fast (~ 390 km s _1 ) into a fully 
ionized interstellar medium. Its proper motion is perfectly aligned with its very elongated 
X-ray nebula, pointing to a physical link between the nebula and the space velocity of the 
neutron star. However, the standard scenario of a bow-shock synchrotron-emitting nebula 
does not explain the phenomenology of our target. Deep radio polarization images could 
help explain the nature of this PWN, by tracing the underlying large-scale magnetic field 
configuration. Unveiling a population of energetic, nearby pulsars missed by radio surveys, 
the Fermi mission is opening new avenues to detect unusual PWNe which could point to 
different ways to deposit pulsar rotational energy into the interstellar medium. 

RPM thanks Andrew Gosling for his support during the preparation of the Gemini 
observations. The research leading to these results has received funding from the European 
Commission Seventh Framework Programme (FP7/2007-2013) under grant agreement n. 
267251. This research used the facilities of the Canadian Astronomy Data Centre operated 
by the National Research Council of Canada with the support of the Canadian Space 
Agency. 
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Fig. 1. — Results of relative astrometry. Two independent, first-epoch images (2009 October 
25 and 2009 October 26) are compared to a second-epoch image (2011 December 24) in the 
upper and middle panels, respectively. Coordinate residuals after image superposition are 
shown for reference sources (blue dots) as well as for the pulsar counterpart (red star), in units 
of statistical errors. The displacement of the pulsar is apparent. The bottom box displays 
the comparison of the two first-epoch images showing, as expected, no displacements for any 
source. 
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Fig. 2. — The field of PSR J0357+3205 as seen by Chandra in the 0.3-8 keV energy range. 
Data from the two observations performed in 2009 have been used (77 ks exposure time). 
The image has been rebinned to a pixel size of 8" in order to ease visibility of faint diffuse 
structures. The large, elongated nebula is apparent. The pulsar counterpart is marked, 
together with the direction of the proper motion (arrow) and the la errors (dashed lines). 
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Fig. 3.— Inner portion of the field of PSR J0357+3205 as seen by the Gemini/GMOS 
instrument in the H a band. The position of the pulsar is marked by a circle (Of! 6 radius). 
No diffuse emission is seen, related to the proper motion direction of the pulsar, marked by 
an arrow. The upper limit to an undetected bow-shock nebula is 5 x ICT 18 erg cm -2 s _1 
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